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Abstract 

(N ■ 

The electromagnetic fields generated by a ring current around a Kerr black hole 
have been found. The acceleration of a charged particle by a force electric field 
along the rotation axis is investigated in the constructed model, as applied to the 

astrophysics of quasars. 

1 INTRODUCTION 

oo 1 

Studying the interaction of electromagnetic fields with the gravitational field of 
a rotating black hole is of great importance in understanding the astrophysics of 
| quasars. Quasars manifest themselves as compact, intense sources of electromagnetic 
radiation, which occasionally have huge narrowly directed jets and which are most 
likely active galactic nuclei. The BlandfordZnajek process [1] is one of the models 
that accounts for observational manifestations of quasars. A magnetohydrodynamic 
£jT)| (MHD) model of plasma accreting onto a rotating black hole underlies this process. 

Through the BardeenPetterson process [2], accretion can proceed only from the 
equatorial plane; therefore, it makes sense to model accretion as a superposition 
of equatorial ring currents. Such modeling is proper if the pair production by an 
induced electric field gives rise to currents that are much weaker than the source 
ring current. The Hawking effect of particle production on the horizon is negligible 
in this model, because it gives a negligible correction when the Compton particle 
wavelength is much smaller than the radius of space curvature. For electrons, this 
corresponds to a black hole with a mass larger than about 10~ 16 Mq. 

In all the cited studies, calculations were performed by using 3 + 1 formalism. In 
contrast to these studies, we use the general covariant formalism of general relativity 
and do not use the approximation of MHD magnetic field line freezing-in in plasma, 
which leads to the condition for the scalar product of the electric and magnetic fields 
being equal to zero (force- free field). We introduce same notation for electromagnetic 
and gravitational quantities as in [3], (§90). [] 

1 Below, we use the system of units =1 — the speed of light, and G = 1 — the gravitational 
constant, for convenience. 
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2 SPECIFYING BOUNDARY CONDITIONS 

Let us specify the boundary conditions that must be imposed on the electromagnetic- 
field tensor components required to determine the latter. First, the classical bound- 
ary conditions must be satisfied: all field components must become zero at infinity. 
Second, the boundary conditions on the horizon must also be added to the classical 
ones. They are required so that, in the frame of reference associated with a freely 
falling observer (FFO), no anomalies arise in its motion as it flies up to the horizon 
due to the electromagnetic-field acceleration (anomalies result in the violation of 
the condition for the FFO electric charge being a test one. Q Therefore, to find 
the conditions imposed on the field components on the horizon, we must write out 
the 4-vector of FFO acceleration and establish which of its components have the 
anomalies on the horizon related to the electromagnetic field. Below, we use the 
coordinates that are at rest with respect to an infinitely remote, static observer. 
The FFO 4-acceleration components are f\ 

W* = w l + (—)F lj u. (1) 



m 



Here, e and m are the FFO test charge and mass, respectively; w l = —Y l km u u'\ 
where T km are the Christoffel symbols; F %1 are the contravariant tensor components 
of the electromagnetic field; and u are the FFO 4-velocity components. Let us 
calculate the 4-acceleration for a radially falling FFO in the Schwarzschild metric. 
The components of the Christoffel symbols and the FFO 4-velocity for this metric 
are given in [3]. Denoting the Schwarzschild radius by r g = 2M, where M is the 
blackhole mass, we have 
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u =-g rr ] u=-y/r B /r; w° = -~{rjrf' 2 ; w 1 = ^(2 + g rr ). 

The remaining u and w l are zero. When approaching the horizon (r — > r ), the 
metric tensor component g rr — > — oo. Therefore, let us write out the asymptotics of 
the squares of the observed 4-acceleration components in the principal (in g rr ) order 
expressed in terms of the contravariant field components: 



W W° -> -g rr (£f r f (F r0 ) 2 ; W r W - g rr g rr + (£) F 
WeW° - gl (^) 2r f 9e0 (F^f ; W V W - <£ (^f r fg vv (F< 



(2) 



We thus see that the singularity in the radial acceleration component does not result 
in the violation of the condition for the charge being a test one for the following 
reasons: (1) the gravitation near the horizon in the radial direction acts a factor 
of g rr more strongly than the electromagnetic field; (2) the zero W component has 
a weak singularity; and (3) the tangential acceleration components have a strong 

2 A more detailed discussion of this requirement can be found in [4, 5]. 

3 Unless stated otherwise, the Latin indices run the series 0, 1, 2, 3, while the Greek indices run 
the series 1, 2, 3. 
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singularity, which can significantly change the FFO trajectory near the horizon 
and can violate the condition for the FFO charge being a test one: e <C ra. The 
contravariant components of the tangential magnetic field must be set equal to zero 
on the horizon, lest this happen. Similarly, the asymptotics of the squares of the 
4-acceleration components dependent on the covariant field components is 



W W° -> -g rr (i) 2 2l» (F r0 ) 2 ; " W r W r -> g rr g rr + (£) F, 



WW* - £ fe) (« 2 ; W^* -> £ fe) <T (F t 



rO 

2 



(3) 
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We thus see that the covariant components of the tangential electric field with the 
strongest singularity must become equal to zero on the horizon. However, the Kerr 
field rather than the Schwarzschild field is of physical interest. Let us write out the 
Kerr metric and its determinant in a nonrotating (relative to remote stars) frame: 

ds 2 = (1 - r -f )dt 2 -£dr 2 - p 2 d6 2 - 

- (r 2 + a 2 + r -*jj£- sin 2 9) sin 2 9 dip 2 + ^ sin 2 edipdt, (4) 
—g = p 4 sin 2 9 

(here, p 2 = r 2 + a 2 cos 2 9, a = r 2 + a 2 — r g r, and a is the Kerr parameter). In this case, 
in the Kerr field, both the space itself and the FFO are drawn into rotation when 
approaching a black hole. The conditions for the tangential electromagnetic field 
components in the Schwarzschild field then change to the conditions for the same 
components in the frame of reference observers in the Kerr field. Expressions (2) 
and (3) give a summary of the boundary conditions for the electromagnetic field on 
the horizon in the Kerr field: 

F" eo ^0; F% ^0; F" rd ^0; F" T<P ^ . (5) 

Let us write out formulas that can be of use in the subsequent analysis. Designating 
K = -9 Jg vv , K + Vl 2 = a sin 2 0/^, we then have 

g =l/g M , g =l/g„,g = - , (] 

y Asin 2 e ' y A sin 2 " 

3 PASSAGE TO THE REFERENCE FRAME 

Since the current is axisymmetric and stationary, the system has two Killing vectors: 
along the translations in time and along the translations in angle (p. Therefore, all 
fields in frame (4) are stationary. In general, this stationarity is not conserved when 
passing to a frame of reference rotating with an arbitrary angular velocity. To show 
this, let us write out the coordinate transformations to a frame of reference rotat- 
ing with a spatially nonuniform angular velocity u(r, 9) required in the subsequent 
analysis: 

g = dx ' k {SI + uSffi + tfySldfv] (7) 

4 In what follows, by the reference frame we mean a frame in which g Q — (see [4, 5]). This is 
the frame in which the reference observers rotate with angular velocity Cl = —g^/g^^ relative to 
remote stars and have a zero angular momentum. 
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The corresponding transformations of the contravariant field components (see [3]) 
are 

piaO __ paO.^ n j?IO(p _ pOcp _|_ fp°^Q aUJ ^ J?' a <P — ]7 a <P _ ujF a0 ; F' rd = F rB '. (8) 

For the covariant components, we derive Q 

F'oup = F atp \ F' aQ = F a0 + ujF aip ] Fq^ = Fq v ; F' re = F r $ + t e a ^ v F atp d^uj . (9) 

In expressions (8) and (9), a runs the values of r and 9. We see from these expressions 
that only the frames of reference rigidly rotating relative to another stationary frame 
(relative to remote stars) are stationary. The stationary frame of reference rotating 
with the horizon angular velocity Q H coincides with the FFO frame on the horizon. 
Therefore, the boundary conditions in this horizon frame, which below is marked by 
a tilde, are the same as (5): 

Feo^O; F v0 ^0; F r0 ^O; F rip - . (10) 



4 MAXWELL EQUATIONS 

Suppose that the current has a density with a deltashaped distribution function in 
the meridional plane: 

f(r, 6) = J% [S(r - r )5(9 - 6 )] /V=g. (11) 

Since the other currents are assumed to be negligible and since the frame is station- 
ary, the toroidal electromagnetic- field components are zero. Therefore, the axisym- 
metric Maxwell equations for the covariant field components outside the horizon 
are 



e a ^d p F a(p = ; 



(12) 



We thus see that the covariant tensor components of the electromagnetic field can 
be represented as 

F a!fi = d a A v ; F a0 = d a A . (13) 

Here, Ai — are the covariant components of the 4- vector electromagnetic-field com- 
ponents. Let us now write the axial Maxwell equations for the contravariant field 
components without electric charges: 



1 d a (^F a0 ) = 



The right-hand part of the second equation in (14) is given by expression (11). 
5 e Q/37 _ 6q ^ — i s the LeviCivita symbol. 

6 A tilde denotes the electromagnetic field components in a stationary frame of reference that 
passes to the horizon reference frame, i.e., to the frame rotating with the FFO angular velocity on 
the horizon. Q„ = „ g „ . 
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5 DETERMINING THE MAGNETIC-FIELD COM- 
PONENTS 

Let us determine the magnetic field in the Schwarzschild metric (a = in [4]). 
Because of axial symmetry, the vector potential of the magnetic field has the toroidal 
component alone. According to (13), we have for its covariant part 

^ = <?VH' ( 15 ) 

Denoting the tensor components in Euclidean space by the subscript o, we then 
have, according to (6), 



F 7V = gfiS} - ^SfflWdfA, = F^{5} - ^Sffi). (16) 



We see from (15) and (16) that if the function Aj is smooth and has no singularities, 
then the boundary conditions (5) are satisfied with the required asymptotics (2). 
The second equation in (14) can then be rewritten as 

L A(V=^T) = ^3 tot- (17) 



where we designated jf ot = j v + ^s9 r (rF^ v ). Next, let us introduce the physical 
vector components by the definition (see [3]) 



B a = B^\g a/3 \ (18) 

In these components, Eq. (17) has the form of the Poisson equation in Euclidean 
space for the vector potential A 1 with source j^ ot . The solution of this equation is 
known to be (see [3]) 

~ 7 A 

^ tlP 00 7V TV 

A a (rfi) = 8«e* / J J ot 2 r' 2 dr' | sin 6' d6' | dip' | . (19) 

Here, e 7 - - is a unit vector [] in the direction of angle <p, \f — r'\ 2 = r 2 + r' 2 — 
2rr'(cos 6 cos 9' + sin 9 sin 9' cos <p'), 9 — is the inclination of vector r to the z axis, 9' 
is the inclination of vector r' to the z axis, and <p' is the angle between the projections 
of vectors r and r' onto the plane perpendicular to the z axis. The integration 
is performed in Euclidean space outside the sphere of radius r g . Solution (19) is 
obtained by iterations; we assume that r g = in the initial iteration. 
Let us determine the magnetic vector 

H\ = -lElexfrF^, (20a) 



or 



F al3 = ^H x e Xal3 . (206) 



J tot e i = 3 tot r sm d cos V 
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According to (18) and (20), the components A v , H r , and Hg are given by 

1 



r sin 6 



^d^rA"). (21) 



For the initial iteration, according to (21), the following expressions can be derived 
for the magnetic- field components from (19) and because of the presence of delta 
functions in (11): 



H 



(r,0) 



r 



H°e(r,e)/r 
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( cos 8 cos 2 ip'(l+x 2 — x sin 8 cos ip') \ j / 
\ x sin 8(l+x 2 — 2x sin 9 cos ip') 3 / 2 J r 

( x sin 9 cos 2 ip'~ cos <p' \ j I 
\ x(l+x 2 — 2x sin 8 cos ip') 3 / 2 J • 



[22) 



Here, we designated x = r/r . In particular, we obtain 



2vr cos 9 



lim (H° e /r) = 0, lim (H° d /r) = , lim H° r - — : . 

e^o K 61 ' ' *-^ v 61 ' r x 3 ' e-o r r (l + x 2 ) 3 / 2 



(23) 



6 DETERMINING THE ELECTRIC POTENTIAL 



To calculate the electric potential A , we use the formula 

F" = {F pa g'"+F fkp g*'}g* 



(24) 



However, the components F a0 must be expressed in terms of Fp and F^, because 
we know the boundary conditions only for them in stationary frames of reference. 
For the poloidal electromagnetic- field components, using formulas (6), (8), and (9), 
and the expression g m g yi = g at * g y0 = S", Q we obtain 



^a0 



J f 



00 a/3 ~ 



g g "A , p + F^(n H -n)/(K + ^ 



(25) 



where L f = [1 - (Q - SI h ) 2 /{k + tt 2 )}~ has the meaning; of FFO Lorentz factor 
relative to the stationary frame of the horizon. Hence, in the second approximation 
in a* = a/M = L/M 2 , the dimensionless black- hole angular momentum (0 < a* < 
1), the first equation in (14) can be reduced to 



AA = -47T (f {Ha) + U (Aq 



(26) 



Here, A — is the Laplace operator in Euclidean space and H a — are the Schwarzschild 
magnetic-field components derived in the preceding section: 



U, 



(A)) 



j^d r {r 2 d r A Q ), 



f(H a ) 



4ttt 2 
a 



x f (1 + r -f + ^) [2 cos 6H r - Anj^r 2 sin 2 6} - sin 6(2 + r f)E e jr 
The expression in round brackets on the right-hand side of Eq. (26) is an analog of 
8 For a running the values of r and 9, we have g$ a = and g 0a = 0. 
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the electric charge density in the Poisson equation. However, because of the second 
term, the A — dependent function U, it can be solved by the iteration method 
assuming that U = in the initial iteration. According to the boundary conditions 
(10), solving Eq. (26) is equivalent to calculating the potential A , produced by the 
density of electric charge p = (/ + U) around a conductive sphere of radius r g in 
Euclidean space. This external problem for the Poisson equation can be solved by 
the image method [6]. The solution that satisfies the boundary conditions (10) is 

/; T It' \ 00 it n 

P(r>,e>) - - f— ^ r' 2 dr' \ sin 6' d6' | d V ' \ . (27) 

Mr — r\ \r — r \r g /r )|J r g o 

The integration is performed in Euclidean space outside the sphere of radius r g [see 
(19)]. According to (7), we obtain the potential A from (27) 

A = A Q Vt H A ip ; Ap = A^ . (28) 

We can factor aJ v /r outside integral (27). Denoting the magnetic-field strength at 
the center of the system in the absence of a black hole by H = 2irJ v /r [see (23)], 
we derive for A and A : 

A = Int (y, #J • a* • M • F /(2tt) ; A = Int (y, 9^ ■ a* ■ M ■ H /(2ir). (29) 

Here, Int y-^-, and Int (^-, are dimensionless functions, which can be numer- 
ically calculated using expressions (27) and (28). The corresponding results are 
shown in the figure. Some important characteristics of the solution can also be 
established analytically. 

1) . Since the system has a mirror symmetry relative to the inversion in the 
equatorial plane, f\ the electric-field direction on the axis depends on the coincidence 
of the directions of black-hole angular momentum and ring-current moment. 

2) . We see from (27) that A — > when r —>■ r H (the expression in curly braces 
becomes zero). This important result has a simple explanation: 

Since the total charge under the horizon is zero, the electric field must pierce the hori- 
zon in different directions. Therefore, there must be an equipotential surface that 
separates these directions, pierces the horizon, and goes to infinity. Since the poten- 
tial A changes neither on this surface, nor on the horizon (in view of the boundary 
conditions), the horizon potential is equal to the potential at infinity (zero). This 
conclusion has an important implication: the existence of a local extremum for the 
potential A . Similar reasoning also applies to A . It might seem that a nonzero 
electric charge density must arise near the extremum. In reality, however, this is 
not the case: the charge density is determined by the contravariant electric-field 
components [F a0 rather than F a o), while the derivation of (27) is based on the first 
equation in (14), a zero charge density in the entire space. We emphasize that there 
is no extremum in the well-known solution of Wald [7] for a black hole placed in a 

9 This can be seen from (26) and from the expression for / (sin# and Hg are even, while cos 6 
and H r are odd relative to the mirror inversion). 
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. 1: Figure 1. The dependence Int(r/r g} 8) for r /r = 2. 



uniform magnetic field aligned along the symmetry axis of the black hole. This is 
because the field does not vanish at infinity in this solution, and this nonphysical 
boundary condition wipes out the extremum. 



7 THE KINETIC ENERGY OF A CHARGED 
PARTICLE EMERGING ALONG THE AXIS 

The conserved massenergy of a charged particle with a zero angular momentum 
component along the z axis is given by [4, 5] 



E = m o7i 



p 2 A 



\ (r 2 + a 2 ) 2 — a 2 a sin 2 9 



- QA = m 7_ 



oo 
L 1 



(30) 



where r y L — is the Lorentz factor of the particle, and q — is its charge. Consider 
the acceleration of an emerging charged particle along the axis by an electric field. 
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Since qA = —\qA \ in this case, the kinetic energy of this particle at infinity with 
j L , g 00 , and A specified at any point of the z axis is 



E k = E-m = m ( v /^~y i - 1) + \qA, 



(31) 



Since the electromagnetic energy of a charged particle for quasars is much larger 
than the corresponding rest energy, it is convenient to represent the result for A 
as the energy (in electronvolts) that an elementary charge emerging along the axis 
acquires. For magnetic fields H « 10 4 and mass M ~ 10 8 M^s we obtain from (29) 



We have calculated the electric field generated by a rotating black hole that interacts 
with an external magnetic field. As can be seen from the above discussion, this field 
is quadrupole in nature, being actually an analog of the pulsar electric field (unipolar 
inductor model [8]). The situation on the z axis was not chosen by chance. Because 
of strong magnetic fields, the Larmor radii of a charged particle must be of the order 
of the gravitational radius of the system, and the particle acceleration mechanism 
will be effective only in directions close to the z axis, where the Lorentz force does 
not act. The reader may ask a legitimate question: Will a strong electric field 
produce electronpositron plasma near the horizon and will it destroy the force field 
(E ■ H) 7^ 0? At large gradients in A , this can actually happen, and the problem 
should then be solved in the force-free approximation (see review papers [911] on this 
subject). However, as numerical estimates show, the electric-field strength in the 
model does not exceed 10 7 V cm' 1 while a strength of the order of ~ 10 9 V cm' 1 is 
required for the particle production. In addition, in the force-free approximation, 
particles cannot be accelerated along the z axis, because there is no electric field 
in this direction. Besides, in any direction in a forcefree field, the electric and 
magnetic fields equally act on the particle. Therefore, even if its trajectory goes 
to infinity, the particle loses most of its energy in the process. Recall that jets are 
observed in quasars precisely in the directions along and opposite to the axis. For 
this reason, solving the problem under consideration in the force-field approximation 
is of considerable importance. 
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